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Abstract 

A one-dimensional finite element model of a sandwich panel with insert is derived using the 
approach used in the Thomsen model. The one-dimensional model produces results that are 
close to those of a two-dimensional axisysmmetric model. Both models assume that the core is 
homogeneous. Our results indicate that the one-dimensional model may be well suited for small 
deformations of sandwich specimens with foam cores. 



1 Introduction 

The numerical simulation of complicated sandwich structures containing inserts can be 
computationally expensive, particularly when a statistical analysis of the effect of variable input 
parameters is the goal. Simplified theories of sandwich structures provide a means of assessing the 
adequacy of the particular statistical technique that is of interest. 

Theories of sandwich structures can be broadly classified into the following types: 



• First-order theories (see for example, [fl]]). 

• Higher-order linear theories that do not account for thickness change (see for example [2] and 
references therein). 

• Geometrically-exact single-layer nonlinear theories that do not account for thickness change 
(see for example, 0). 

• Higher-order linear single-layer theories that account for thickness change (see for example, 

BH30). 

• Higher-order linear multi-layer theories that account for thickness change (see for example, 

ClIHlEllTDlini) • 

• Higher-order nonlinear single-layer theories that account for thickness change (see for 
example, D2 EH 05 HUES!). 

'Corresponding author, email: b.banerjee@irl.cri.nz 
t email: bryan.smith@irl.cri.nz 



1 



2 THE THOMSEN MODEL 



Most theories start with ad-hoc assumptions about the displacement or stress field. 
Geometrically-exact theories avoid such assumptions but are hampered by the requirement that 
special constitutive models have to be designed for consistency. The linear theory proposed by 
Thomsen and co-workers HO provides a formulation that is simple enough to be evaluated 
rapidly. Therefore, we have chosen that formulation and applied it to an axisymmetric sandwich 
panel in this work. The work of Thomsen involves the solution of a system of first order ordinary 
differential equations using a multi-segment numerical method, We have instead chosen to use the 
considerably simpler finite element method to discretize and solve the system of equations. 



2 The Thomsen Model 

Since we are considering a simplified axisymmetric form of the sandwich panel problem, we start 
with the governing equations expessed in cylindrical coordinates. The geometry of the sandwich 
structure under consideration is shown in Figure [T] 
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Figure 1 - The geometry of the sandwich panel. 



2.0.1 Strain-displacement 

The strain-displacement relations are given by 
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Axisymmetry implies that the displacement ug = ug[r) and all derivatives with respect to 6 are zero. 
If in addition, the displacements are small such that ug = C r (this assumption is not strictly 
necessary), the strain-displacement relations reduce to 
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2.0.2 Stress-strain 

The stress-strain relations for an orthotropic material are 
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In cylindrical coordinates 







"Cn 


C\2 


Cl3 








" 










Cl2 


C22 


C23 



















Cl3 


C23 


C33 




























C44 










£0z 


<J TZ 
















C55 







&rz 


_a r e_ 




_ 














Cm. 




e r e_ 



From axisymmetry, we therefore have 



O rr = C\\ Err + C\2 £09 + Cl3 E zz 

&ee = C12 £rr + C22 £ee + C23 £zz 

&zz = Cl3 Err + C23 ^96 + C33 &zz 
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2.0.3 Equilibrium 

We assume that there are no inertial or body forces in the sandwich panel. Then the 
three-dimensional equilibrium equations take the form 
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The equilibrium equations in cylindrical coordinates are 
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Because of axisymmetry, all derivatives with respect to 9 are zero and also gq z and a r g are zero, the 
reduced equilibrium equations are 
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2.1 Facesheet equations 

The facesheets are modeled using the Kirchhoff-Love hypothesis, i.e., that transverse normals 
remain straight and normal and that the normals are inextensible. In that case, the displacement field 
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in the plate takes the form: 

u r (r, 9, z) = UQ r (r, 9) - z ; u e (r, 9, z) = u oe (r, 9) - z ; u z (r,9, z) = w (r,9) (10) 

where uq t is the displacement of the midsurface in the r-direction, uqq is the displacement of the 
midsurface in the ^-direction, and wq is the z-direction displacement of the midsurface. 
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We define the stress resultants and stress couples as 

rf f 

Gee dz ; M rr := 
-f J-s J-f 

where the thickness of the plate is 2/. 
2.1 .1 Strain-displacement relations 

From axisymmetry, the strain-displacement relations are (for small rotations, i.e., NOT the von 
Karman strains) 
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Plugging in the displacement functions in the strain-displacement relations gives 
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To simplify the notation, we define 
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to get 



e rr (r, z) = e° r (r) + z e l „(r) ; e ee (r, z) = e° ee {r) + z e\ e {r) 
2.1 .2 Stress-strain relations 

Assuming that the facesheets are transversely isotropic and taking into account the 



strain-displacement relations ( 13 1, the axisymmetric stress-strain relations are 
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If we make the plane stress assumption, a zz = 0, then we have 
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Then the relations between the stress resultants and stress couples and the strains are 
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From the expressions for strain in equations ( 15 ) 
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Therefore, the relations between the stress resultants and stress couples and the strain can be 
expressed in matrix form as 
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where A%j = 2/ Cjj are the extensional stiffnesses of the plate and Dij = 2/ 3 /3 are the bending 
stiffnesses of the plate. 

2.1.3 Equilibrium equations 

The plate equilibrium equations may be derived directly from the three-dimensional equilibrium 
equations. However, it is more informative to derive them from the principle of virtual work 

SU = 5V exX (23) 

where SU is a variation of the internal energy and 5V ext is a variation of the work done by external 
forces. 

The variation in the internal energy is given by 
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where Oq represents the reference surface of the plate. In terms of the definitions in ([14]), 
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The variation in the work done by the external forces is 



5V ext = / [q(r) 5w + s(r) (5u 0r - Zf 5(3) + p(r) 5u oe ] d£l 
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where q(r) = g To P Face ( r ) + gBotFace^ r -j j s a distributed surface force (per unit area) acting the 
positive z direction, p(r) = p To P Face ( r ) _|_pBotFacej- r ^ j s a distributed surface force (per unit area) 
acting the positive r direction, s(r) = s To P Face ( r ) -|_ s BotFace^ r ^ j s a distributed surface force (per 
unit area) acting the positive 9 direction, zj takes the value +/ at the top of the facesheet and — / at 
the bottom of the facesheet, and t = t r e r + tg eg + t z e z is the surface traction vector. 

A schematic of the loads thare are applied to the facesheet is shown in Figure [2] 




Figure 2 - The loads on a facesheet. 



In terms of resultants over the thickness of the plate 
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Integrating the 5(3 term by parts over the area Q,q gives 
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Then, from the principle of virtual work, we have 
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Invoking the fundamental lemma of the calculus of variations (and keeping in mind that the 
displacement variations and the applied tractions are zero at points on the boundary where 
displacements are specified), we get the governing equations for the axisymmetric plate: 
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Then the boundary conditions are 
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The governing equations are of order 6 in the displacements (uo r , wq) and there are 6 nontrivial 
boundary conditions, (uo r , wq, dwo/dr, N r , Q z , M r ). 
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2.1 .4 Summary of facesheet governing equations 

The governing equations for the plate can then be summarized as follows: 



• Equilibrium equations: 
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• Strain-displacement relations: 
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2.1 .5 Conversion into first-order ODEs 

We would like to convert the governing equations for the axisymmetric plate into ODEs of first order 
for computational purposes. To do that, we note that the stress resultants are related to the 
displacements by 
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Using (EFT), 
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Then the relations between the stress couples and the displacements take the form 
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Plugging in the expression for Mee from (|5"T|) and the expression for the derivative of p ([52]) we have 
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To reduce the order of the equilibrium equation for the stress couples, (41 ), we note that taking the 
derivative of Q r from ([53]) gives us 
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2.1 .6 Summary first-order ODEs for facesheets 

The ODEs governing the facesheets are: 
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2.2 Core equations 

2.2.1 Stress-strain relations 

We assume that the core is transversely isotropic. In that case, the stress-strain relations in the core 
have the form 
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(?ee = C12 £rr + Cn £ee + Ci3 £zz 

Ozz = C13 Err + C13 Eqq + C33 E zz 

&6z = ; cr rz = C55 £ rz ; u r g = 
If we also assume that the core cannot sustain any in-plane stresses, then 



(65) 



Therefore we have 



= = Cn £ rr + C12 Eee + Cl3 &zz 

= = C12 Err + Cn £ee + Cl3 £zz 
(Cn - C12) (^r - eee) = 



(66) 



(67) 



which implies that Cn = C\2- If we assume that Cn = C\i = e C13 where e <C 1 is a positive 
quantity, then we have C13 = 0. Therefore the stress-strain relations in the core reduce to 
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2.2.2 Strain-displacement relations 

From the strain-displacement relations we have 
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2.2.3 Stress-displacement relations 

Using the stress-strain relations we get 
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2.2.4 Equilibrium equations 

The equilibrium equations also reduce accordingly to 
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2.2.5 Expression for u z 
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where A(r) is a function only of r. Integrating by parts, we have 
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Now we assume that the displacement u z is quadratic in z to get 
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At z c = c the displacement of the core is equal to the displacement of the top facesheet, i.e., 
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We can also calculate the displacement at the bottom facesheet 



w 2 {r) = u z (r, -c) = S : 



;>:>, 



-c a zz (r, -c) - B(r) 



+ A(r) 



Again, eliminating A(r), we have 

w 1 ^) - w 2 (r) = c S33 [a zz (r, c) + a zz (r, -c)] . 

2.2.6 Eliminating a zz 
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We would like to eliminate a zz from the expression in equation ( |90| >. To do that, we recall that 
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cr zz (r, z c ) = [w l (r) - w 2 (r)] + z c B(r) 



(81) 
(82) 

(83) 

(84) 
(85) 
(86) 
(87) 

(88) 
(89) 
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2.2 Core equations 



2 THE THOMSEN MODEL 



Using (188]) and ((89) in d78) gives 



u z (r,z c ) = w l (r) + 



Zr. - C 



w 1 (r) — w 2 (r) 



+ (*c + c) S33 S(r) 



Now, from equations ( 10 1 for the facesheets, we have 

w 1 (r) = w t ° p ; w 2 (r)=w b ot 



respectively. Plugging these into (90 1 gives 



u z {r,z c ) = \ ^+1 -r-Hf-l -o bot + ^ « " c 2 ) 5 33 B(r) . 



(90) 



(91) 



(92) 



2.2.7 Expression for u T 

Recall that 

Therefore, 



C55 



3u r <9-u 2 
<9z dr 



du r du z 

= 1 055 o>z — ^— ; 
oz or 

Also, taking the r-derivative of equation d92b, we have 



S55 : = 1/C; 



55 



1 / 2c 

2 c 



+ 1 



gwg 1 / z c 
dr 2 \ c 



dw. 



bot 



dr + 2^~ C ' 



dB 
dr 



(93) 



(94) 



(95) 



Substitution of (95 1 into (94 ) gives 



du r 
dz 



2 S55 cr r - 



1 / £c 

2 c 



+ 1 



+ 2 7" 1 "^"2^- C ) 533 



dr 



d(J r 



Note that ^ z = ==>■ cr r2 = er rz (r). Integrating (96) between and z c , we get 

u r (r,z c ) = 2 S55 z c (T, 



\ dw l Q P 

~2 \2c + Zc ) ^T + 2\2c~ Zc 



1 z 2 c \ dw l ° p 1 (zl 



dw% ot 
dr 



1 / 4 9 \ d# 

2 I- C2 H 5 ^ + G « 



(96) 



(97) 



where G(r) is a function only of r. 
At z c = c,u r = u 1 (r). Hence we have 

G(r) = u 1 - 2 S55 c a rz + 



3c dvjQ V c dw, 



+ 



hot 




4 dr 4 dr 



c 3 dB 
J S33 ~dF 



(98) 
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2.2 Core equations 



2 THE THOMSEN MODEL 



Substitution of d98]> into d97l) gives 



u r (r, z c ) = u 1 + 2 5 55 (z c - c) o rz + 



3c I lz 2 c 
T~ 2 I 2~c 



+ £ c 



dr 



+ 



c l(K 

4 + 2 I 2c 



c° 1 / z, 
3~ + 2 \ 3 



c 2 

C Z c 



c^ 0t 
dr 

dB 



S33 



dr 
(99) 



Now, from equations ( 10 1 and (30) for the facesheets, we have 

dm top dm bot 



11 



dr dr 
1 (r) = C + / t0P P° P ! «V) = «o? - / b0t /3 b0t 



(100) 



where 2 / top and 2 / bot are the thicknesses of the top and bottom facesheets, respectively. Plugging 
these into d99] gives 



u r (r, z c ) = < r p + / top /3 top + 2 S 55 (z c - c) a rz + 



3c 1 / z 2 c 
T~ 2 I 2c 



+ z c 



+ 



c 1 / z£ 
4 + 2 I 2c 



bot 



■ + 



1 Me 



3 2 V 3 



c 2 z c 



^33 



/3 top 

dB 

dr 



(101) 



or, 



u r (r, z c ) = itQ 0p + 



3c 1 zi 



+ z c 



4 2 V 2c 

2 5*55 (z c - c) cr r 



/3 top + 



c 1 / z% 
4 + 2 I 2c 



■ + 



2 I 3 



c 2 z c 



•933 



/3 bot 

dB 

dr 



(102) 



2.2.8 Governing equation for the core 



Now, at the bottom of the core, z c = — c. From ( |102[ ) we have 

2 dB(r) 

u r (r, -c) = u 2 = < p + [/ top + c] /3 top + c /3 bot - 4 S 55 a rz c - - S 33 c 3 (103) 



Also 



u 2 = u b f - / bot /3 bot 
Hence 

2 dB(r) 

= u£P - u b ° l + [/ top + c] /3 top + [/ bot + c] /? bot - 4 S55 a„ c - - S 33 c 3 -A^ (104) 
Recall that 

dd rz cr rz 
dr r 

Therefore, 

(f5 d 2 

dr (ir 2 r 2 r dr 
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3 COUPLED GOVERNING EQUATIONS OF THE FACESHEETS AND THE CORE 



Plugging (|105|> into (|104|) gives 



= < p - u b f + [/ top + c] /? top + [/ bot + c] /3 



hot 



4 555 o> 2 c - - S 33 c 



+ — 



dr 2 



r dr 



or, 



+ 



<ir 2 r dr 



1 6^33^55 \ 
9 H 9 <7rz 



3C 33 r 

' 2c 3 



top hot 

u 0r u 0r 



+ (/ top + c) /? top + (/ bot + c) /3 



hot 



2.2.9 Conversion into first order ODEs 



To convert ( 107 1 into first-order ODEs, we define 



T r :-- 



da r 



dr 



Then equation ( 107 1 can be written as 

3C33 



dT r 
dr 



, top _ hot 
a 0r u 0r 



2c 3 



+ (/ top + c) /? top + (/ bot + c) /3 



hot 



(106) 



(107) 



(108) 



(109) 



2.2.1 Summary of first order ODEs for the core 

The governing equations for the stresses in the core are 

do rz 



dr 
dT r 
dr 



cr Z z(r, z c ) 



T r 



3C: 



2c 3 



11, 



top _ ^bot 
Or 



<+{f* + c) /3 top + (/ bot + c) (3 



jbot 



1 6C 33 S 55 \ 

9 H 5 0>z 

r z c z I 



r 



C33 
~2cT 



\w x (r) - w 2 (r)] - z c T r - z t 



(110) 

(111) 

(112) 



3 Coupled governing equations of the facesheets and the core 

In the previous section, ODEs have been derived that partially couple the core to the facesheets. To 
complete the coupling of the facesheets to the core we have to balance the forces at the interfaces 
between the core and the facesheets. We introduce some new notation to aid us in the coupling 
process. Recall that for a facesheet 



s ^ — ^TopFace _|_ g BotFace . q(j-} _ ^ Top Face _j_ ^BotFace 



(113) 
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3 COUPLED GOVERNING EQUATIONS OF THE FACESHEETS AND THE CORE 



We identify these two sets of applied tractions on the two facesheets using the notation 

s to P(r) = s" + s tb ; c/ top (r) = q tx + q tb ; s bot (r) = s bt + s bb ; q bot (r) = q bt + q bb 

The tractions at the core-facesheet interface are given by 

t = t r e r + tg eg + t z e z 
= (n r a rr + ng a r g + n z a rz ) e r + (n r a r g + ng agg + n z ag z ) eg 
+ (n r arz + ng ag z + n z a zz ) e z 

where e r , eg, e 2 are the basis vectors in the r, 6, z directions. In the core 
a rr = agg = og z = a r g = 0. Therefore, the traction vector simplifies to 



(114) 



(115) 



t = n z a rz e r + (n r a rz + n z o zz ) e z (1 16) 

At the interface between the core and the top facesheet, n r = 0, n z = 1 while at the interface 
between the core and the bottom facesheet n r = 0, n z = — 1. Therefore, 

t tc (r) = < re (r) e r + <x™ re (r, c) e z ; t bc (r) = -a™%r) e r - a c z T(r, -c) e z (117) 

To couple the facesheet equations to the core equations we have, due to the continuity of tractions at 
the core-facesheet interfaces, 



tb 


V) + t te (r) 


• e r 


= = 


=> s tb {r) = 




bt 


V) + t bc (r) 




= 


s b \r) = 




tb 


V) + t tc (r) 


■ e z 


= = 


* q tb (r) = 




bt 


;r)+t bc (r) 


■ e z 


= 


q bt (r) = 


-c) 



(118) 



From equations ( |89| ) and (|7f 

C(r,vl 



C33 
~2c~ 



top 



W Q 



bot 



da* 



dr 



(119) 



Therefore, 



~tb/ 



,bt/ 



a 



33 



2c 



if, 



top 



?/i bot 



""rz u rz 
+ C : h C 



33 



2c 



■»', 



top 



..bot 



+ C 



dr 



dr 



+ c 



CT, 



r 

core 



Equation (57 ) then takes the form 



(120) 



dQT _ 




q top (r) = 




dr 


r 


r 


dQ bot _ 




q bo \r) = 


Q b r 0t 


dr 


r 


r 



+ 



a 



33 



2c 

/^(core 
°33 

2c 



top 



bot 



top 



w 



bot 



da c r °J e 

dr 
da c r f e 

dr 



Similarly, equation ([49]) takes the form 



dN, 



top 



fir 



.top .top' 

^12 ^11 



4? 



top 



+ 



(4T) 2 - 



top 

Or 



r 

T core 



,bb 



2 ^rz 



core ^tt 



(121) 



(122) 
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3 COUPLED GOVERNING EQUATIONS OF THE FACESHEETS AND THE CORE 



and 



dN 



bot 



dr 



A bot A bo 
^12 ^11 



N, 



hot 



Also, equation (55 ) takes the form 



dM, 



top 



dr 



n top n top' 
-^12 -^11 



and 



dr + 



D top 

^11 



nbot nbot 

-^12 ~~ -^11 



top 



■ + 



■ + 



A\f 

PIT) 2 - PIT) 2 



bot 

Or core _ bb 

„2 * 



^11 



^top 



nbot 
^11 



bot 



Pi?) 2 ~ (Ujf ) 2 



r)bot 
^11 



/3 ! 



bot 



(123) 



(124) 



.rbot /^bb , ^.core\ /-ioc\ 
2 / I s +^2 ) ( 125 ) 



The governing first order ODEs for the facesheets and the core can then be expressed as 



Top facesheet: 



du 0r 

dr 

dw l Q P 
dr 

d/3 to v 
dr 

dN* p 
dr 

dM?? 
dr 

dQr P 
dr 



N top .top top 
^top ~~ ^top r 



pop 

top 
12 



D$° p /3 to P 



.top .top 
^12 ^11 



47 

Q? p + 



■ + 



«iT - (4°l) 2 



A*? 



top 

x 0r 



core tt 

a vy - S 



n top n top' 

-^12 -^11 



top 



D 1 



Q to P /'-'core 

r 2c 



M, 



top 



+ 



PIT) 2 - Pi°f 



top\2' 



-,top 
[1 



< P (r) - O) 



D 1 



+ / to P(s tt_ a cor e) 



(126) 
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4 FINITE ELEMENT FORMULATION OF THE COUPLED GOVERNING EQUATIONS 



Bottom facesheet: 



Core: 



i, bot 
""Or 


Arbot 
rr 


A bot „ bot 
^■12 "Or 


dr 




4 bot 
^11 


4 bot 
^11 


r 




rlin hot 

LiU/Q 


= /3 bot 
M b r ot 








dr 
d(3 hoi 


^bot ^bot 


dr 




rjbot 
^11 


rjbot 
^li 


r 






" A bot 


4 bot 


TW-bot 


dr 




4 bot 
^11 




r 


dM b r ot 


= Q** + 


n bot 

^12 




rjbot 
-^11 


dr 


nbot 


dQ* ot 




Q b r° l 


/'"'core 
°33 


top 

w 


dr 




r 


2c 





+ 



(ASf)a _ (A bot )2 



M 



bot 



+ 



A\° 
(D 



,bot 

i 0r 



„bb 



bot\2 
12 J 



bot\2 
11 J 



D bot 







bot 



j?bot / „bb i ^co 



a,. 



-,bb 



(127) 



dr 
ciT r core 
dr 



rpcore 



2£<core 



top 
x 0r 



t^+C/^ + c) /3 top + (/ bot + c) /3 

rricore 



bot 



(128) 



+ 



g^core^core 



This is a set of 14 coupled ODEs that can be solved using a number of approaches. Thomsen and 
coworkers QUI use a multi-segment integration approach to solve these equations. Since it is 
considerably simple to solve the original system of equations using the finite element approach, we 
have used finite elements in this work. 



4 Finite element formulation of the coupled governing equa- 
tions 

For the finite element formulation of the governing equations, it is convenient to start with the 
statement of virtual work for the facesheets, i.e., 



where 



a, 



d5u t 



Or 



dr 



H SU() r — M rr 

r 



d 2 5wr i Mqq d5wQ 



dr 2 



r dr 



dQr 



n 



d5wo 

s(r) 5u 0r - z f s(r) + q(r) 5w 



+ 



dSwo 

N r 5uQ r — M r — h Q z 8wo 

dr 



d£2r 



dTr 



du, 



Or 



U 0r 



du, 



Or 



f A 



u 0r 



ii 



N rr = A u —— + A 12 — ; N ee = A 12 — 

dr r dr r 

d 2 Wn D\o dwn d 2 wn Du dwn 

M rr = -D ll —^- — —^; M ee = -D 12 



dr 2 r dr 



dr 2 r dr 



(129) 



(130) 
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4 FINITE ELEMENT FORMULATION OF THE COUPLED GOVERNING EQUATIONS 



Separating terms containing 5uo r and 5wq leads to two equations 



Qo 



d5u 0r 

N rr + 



d 2 5w \Mqq 



dr 



■ + z f s{r) 



s(r) > Suor 



ddwo 
dr 



+ q(r) 8wq 



d£lo 
dflo 



N r Suor dTo 



r 



r 



d5wo 

M r — Q z 5w 

dr 



(131) 
(132) 



The continuity of tractions across the facesheet-core interfaces requires that 



sjcore 

q top (r) = q tb (r) + q*(r) = -- [w ° p (r) - w h ot (r)} + c 



2c 

rjcore 

7 bot (r) = q U (r) + q bb (r) = -f- [w^(r) - w bot (r)} + c 



rz , v rz . tt/ 

^ + c +q n rj 

JfT core (T COre 

^ + c ^ + (? bb (r) 

ar r 



(133) 



and 



s to P(r) = s tb (r) + s tt (r) = -a™ re (r) + s» 
s bot (r) = s bt (r) + s bb (r) = <° re (r) + s bb (r) 
Plugging these into equations ( 1 3 1 [ > and (132]) leads to, for the top facesheet, 



(134) 



N top 



n 



p dSu 0r 


+ 1 




dr 


r 


d 2 5w Q op 






dr 2 


l r 




-1 


C /^(core 

[* < 



(7 



core _ „tt I fatop 



dQ = <b iV r top 5u Q op dT (135) 



/t op ( s tt 



(i5w;Q 0p 
dr 



dr r 

^op^T_ Q top ^top 



dr c 



(136) 



and for the bottom facesheet 



bot 



A/ 



hot ^ n 0r 



rr dr 



n 



M bot 



d 2 5 



bot 



dr 2 



■ + 



/V bot 
r 

MS? 



core _ bb 



ou 0r 



i j^bot / bb , ^coreN 
+ / [ s + cr rz ) 



dQo 



dSw^ 1 



N bot 6v$ dT 



v 



dr 



+ 



(~tcore 

(w ° p ~ ^ 0t ) + c 



rz " rz 
+ C 



rir 



M bot^0__ Q bot ^bot 



dTo 



dVto 



(137) 



(138) 



The governing ordinary differential equation for the core is 



d 2 a c r ° z re 1 da c rz K 



dr 2 



dr 



■ + 



g^core^core x 



a,. 



o/^core 
6L/ 33 

2c 3 



top 
l 0r 



^ t + (/ t ° p + c) 



top 



dr 



°-+(f bot + c) 



dw. 



bot 



(139) 



dr 
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4 FINITE ELEMENT FORMULATION OF THE COUPLED GOVERNING EQUATIONS 



Multiplying the equation with a test function and integration over the area yields, after an 
integration by parts, the equation: 



dcr™ re d5a r 



dr dr 



+ 



1 



dSa„ da c r f e 

H ; — 8a T 



dr 



dr 



+ 



i 6qg re s^ re 



a,. 



6U 33 

2c 3 



dm top Hm bot 

C P -*+(/ top + c) ^+(/ b0t + c) 



5a r2 



dtlr 



rz 



dr 



5(T TZ dTo 



(140) 



Equations ( 135| >, ( 136| >, ( |137 1, ( 138 ), and ( 140) form the system that has been discretized using the 
finite element approach. 

We assume that the fields u ^F, Ug° l , u>q° p , w^ 1 , <r™ re can be expressed as 

nu nu 

\r) = < P N?(r) ; u^(r) = £ u\ ot iV«(r) 
i=i i=i 

nw nw 

i=l 
ns 



i=l 



(141) 



i=l 



where nu, nw, ns are the number of nodes and N™' w ' s are the basis functions that are required to 
represent the field variables. Then, the stress and stress couple resultants can be expressed as 



i=l 



bot 



N, 



Ml ? 



M, 



bot 



E 

3=1 



Ibot 
11 



dr 


+ <l 


Nf 
r 


dNf 
dr 


+ <2 t 


Nf 
r 



top , ATtOp 

V /)/) 



?! 



bot . Arbot 



E 

nu 

E 



.top » JV J . .top iv j 

iZ dr 11 r 



.4 



, t dNV , . Nf 

bot J , a bot J 

^ dr 11 r 



top 



hot 



E 

i=i 

nw 

E 



n to P d 2 ^7 -Djf dA 7 

-^11 . n + 



dr 2 



top 



. t d 2 Nf D h ° l dNf 
nbot j . iz j 

u ll r-s- + 



dr 2 



dr 



,„bot . ii/rbot 



E 

i=i 

nw 

E 



J2 ATM rjtop ,nr W 
top U iV J -k'll " iV J 



12 



dr 2 



dr 



w 



top 



j2aj«i r)bot A M w ' 
bot " JV J , ^11 " iV j 



12 



dr 2 



■ + 



dr 



(142) 
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4. 1 mitFMWEnEbMMEMTtmRMULATION OF THE COUPLED GOVERNING EQUATIONS 



and the momentum balance equations can be written as 



/ 

/ 

/ 



N top <^\_ . 
rr dr ^ 

rr dr 

M to P 
rr dr 2 



fc=l 



dflr 



3 tt iV^ dn + / ^V' 0P N i dT (143) 



N. 



bot ns 



jr a fc ivA iV^ dft = / s bb iV? dft + / iV r bot N? dV (144 



M, 



top 



+ H^ top I> iv * 



fc=i 



dm 



fc=i 



dr r 



= / / top S tt -^^o + 
•/n dr ^r 



n 



d 2 N w 



dr 
M bot 



dm 
M to P «^ 



dr 



/^core 
°33 

2c 



EK op - ^ ot ) ^ 



k=i 



dTn 



dr 2 



f/ bot J2^ N k\ 
k=i ) 



dm 



k=i 



dN s k Nl 
dr r 



-J 

Jn 

r r ns ( 
Jn 



d7V w 

/ bot a bb^ ! _ d j 2() + 
dr 



dr 
df^o 



M bot - Q bot AT™ 



^- ek op - ^ ot ) w 



fc=i 



divi N?\ dm dm m 



dr 



3 _|_ J 



dr 



dr r 



dr 



1 6C|§ re 5§^ re 



dr 



iV| Nf \a 



^33_ 

2c 3 



diY ; 



,fc=i 



/■- _|_ rz 



dr 



fc=i 

iv/dr 



E ^ M « op - ^ ot ) + E (/ top + c ) «4 op + (/ bot + c ) ^ ot 



(145) 



(146) 



d^o 
(147) 



After plugging in the expressions for the resultant stress and stress couples, we can express the 
above equations in matrix form as 





[0] 


[0] 


[0] 


[kL]" 




- u top- 




r f top-i 


[0] 


fK" 1 

L ww J 


[0] 


[0] 


[kL] 




w to P 




flop 
1 tu 


[0] 


[0] 




[0] 


[kL] 




u bot 




fbot 

u 


[0] 


[0] 


[0] 


\K bb 1 






w bot 




rbot 
i w 


JK*J 


[Kw\ 


[kL] 


[kL] 


[K M ] 











(148) 



4.1 Finite element basis functions 

Note that the stiffness matrix is not symmetric. This system of equations is solved using 
COMSOL™using quadratic shape functions for the u-displacement and the cr-stress and cubic 
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Hermite functions for the ^-displacement, i.e., in each element nu = ns = 3, nw = 4, and 

N?(r)=N°(r)= LZ^l 
j. j. r r 

j=i&j (149) 
JVf (r) = 1 - 3 r 2 + 2 r 3 ; N??(r) = (r 2 - n) (r - 2 r 2 + r 3 ) 

JV^(r) = 3 r 2 - 2 r 3 ; iV^(r) = (r 2 - n)(-r 2 + r 3 ) 



4.2 Boundary conditions 

The natural boundary conditions are 

Top f acesheet : i\£ op , Q l ° v , M^ op 
Bottom facesheet : iV bot , Q bot , M, 



(150) 



Core: — ^ + 

dr 



The essential boundary conditions are 



duu 

Top facesheet : uf, op , !d1 op , — -^ ) — 

ar 



Bottom facesheet : u^ 1 , w^ ot , 
Core : a', 



dw\\ ot (151) 



core 

rz 



Note that fixing the u z displacement at the boundary of the core is equivalent to setting the natural 
boundary condition in the core to zero when Wq P = Wq 01 = at the boundary. 

4.2.1 Through-the-thickness insert 

The boundary conditions used for a simply-supported sandwich panel with a through-the-thickness 
insert are shown in Figure[3] The radius of the insert is r^, that of the potting is r p , and that of the 
panel is r a . Therefore, for part of the panel, the potting is assumed to have the same behavior as the 
core. 

For this situation, the boundary conditions at the left edge, r = r,, are 

= u b f = ; M^°p = M r bot = 

^top _ ^ J . rfxA _ ^ J 

Wz /top + 2c + / bot ' z f°P + 2c + / bot ( 152 ) 

Q 



core 

A' 

where A = 2 ir ri (/ top + 2c + / bot ). These are applied to the two dimensional model as a constant 
pressure in the z direction, with P = Q/A. 

At the right edge, r = r a , the structure is simply supported, with the conditions: 

w ^P = w bot = Q . M top = M bot = . ^top = N bot = q ^ + f™ = fJ> ^ 

dr r a 

The support condition is applied to the two dimensional model by setting wq = along the right 
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4.2 murBBtifTErMlEMENT FORMULATION OF THE COUPLED GOVERNING EQUATIONS 
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Figure 3 - Boundary conditions for a through-the -thickness insert in a simply 
supported sandwich panel. 



edge. 
4.2.2 



Potted insert 



The boundary conditions used for a simply-supported sandwich panel with a potted insert are shown 
in Figure[4] The radius of the insert is rj, that of the potting is r p , and that of the panel is r a . The 
length of the insert is 2/, and the thickness of the potting below the insert is 2c~. Therefore, the 
insert is being treated as a thin plate in the region above the potting and the potting is being treated 
as a material with features similar to the core. 

To allow for the jump discontinuities on the two sides of the insert-facesheet interface, we define the 
quantities u^ r and u J, to be the u$ T displacements of the insert and the top facesheet, respectively. 
The locations where these quantities are evaluated are shown in Figure |4| Then the continuity of 
displacements requires that 



''Or 



Or 



(fi ~ f) 



^interface 

dr 



(154) 



There is also a jump in the shear stress in the two sections of the potting to the left and right of the 
interface. Let these quantities be a~ z and cr+,. We assume that the average force at the interface is 
balanced, i.e., 



c a' 



(155) 
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5 MODEL TEST CASES: FRP SANDWICH 




The boundary conditions at r = are 

top hot n . UW UW n . n top n bot n . core n n ^ 

u or = u or = u ; —j^r = = u ; Qz = Qz = u ; ^ = u ( 156 ) 

The simply-supported boundary at r = r a once again requires that 

w top = w hot = ; M' op = M r bot = ; 2V? P = iV bot = ; ^ + ^i = . (157) 

ar r„ 



5 Model Test Cases: FRP Sandwich 

In order to validate the one dimensional approximation, the results for test cases are compared with 
the results generated by a two dimensional axisymmetric model. In each test case, a rigid, through 
the thickness insert applies a vertical compression load of Q = 1000N to a simply supported 
sandwich panel. 
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5.1 Example 1: Stiff facesheets 



5 MODEL TEST CASES: FRP SANDWICH 



The facesheets are assumed to be isotropic, i.e, 



A u = 2 / Cn = 2/ — ^ ; A 12 = 2fC 12 = 2f 



1_„2' ^ "* " J \- V 2 

2 f 2/ 3 E 2/ 3 2f 3 uE 

Dn = ~tt~ Cn = — - 5 ; D12 = -—— C\ 2 



(158) 



3 ~ Xi 3 1-i/ 2 ' iZ 3 iZ 3 1-i/ 2 
where E 1 is the Young's modulus and v is the Poisson's ratio. 
The core moduli are given by 

C 33 = -^; C 55 = 2G h (159) 

where Eh, Gh, Vh are the Young's modulus, shear modulus, and the Poisson's ratio of the core. The 
potting moduli are also computed in a manner similar to those of the core. 

5.1 Example 1 : Stiff facesheets 

The first example problem is taken from |7 ], with the parameters given in Table [T] Figure |5Ja) 
compares the resulting out of plane displacements from the sandwich theory and the two 
dimensional axisymmetric simulations. While there is a small amount of disagreement in the potting 
region, the overall results match up well. The radial displacement, u r , is shown in figure |5Jb), and 
these results match as well. Core shear stresses and transverse stresses at the bottom of the core are 
shown in Figure|6] The stresses match reasonably well too. 

Table 1 - Geometric and material parameters for Example 1 

Geometry (mm): n = 7.0 r p = 10.0 r a = 60.0 c = 10.0 / top = 0.1 / bot = 0.1 

Face Sheets (GPa) : E x = 71.5 G x = 27 .5 

Potting (GPa) : E p = 2.5 G p = 0.93 

Honeycomb (MPa) : E h = 310 Gh = 138 
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Example #1: Transverse Displacement - Top Face 



■ ID Approximation 

- 2D Axisymmetric Model 
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60 
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60 



Figure 5 - Comparisons of displacements from one-dimensional and two- 
dimensional finite element simulations for the model in Table^ 

5.2 Example 2: Soft facesheets 

The second example problem is taken from [8], with the parameters given in Table [2] Once again, 
figure [7J a) shows the out of plane displacements given by the sandwich theory and the axisymmetric 
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6 S UMMARY AND CONCL USIONS 
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Example #1: Core Shear Stress 



Example #1: Transverse Stress - Core Bottom 
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ID Approximation 
2D Axisymmetric Mode: 
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(b) Transverse stress (a zz )- core bottom. 



Figure 6 - Comparisons of stresses from one-dimensional and two-dimensional 
finite element simulations for the model in Table^ 



simulations, and figure |7Jb) gives the radial displacements. As in the first example, the results match 
reasonably well, suggesting that the sandwich theory captures the important physics of the problem. 
The stresses shown in Figure [8] also show that the one- and two-dimensional models predict similar 
results. The values of transverse stress and displacement at the bottom of the core are shown in 
Figure [9] 



Geometry (mm): 
Face Sheets (GPa) : 
Potting (GPa): 
Honeycomb (MPa): 



Table 2 - Geometric and material parameters for Example 2 

a = 150.0 c = 5.0 / top = 0.5 / 



n = 



10.0 
= 40.0 
= 2.5 
= 310 



r p = 

G p 
G h 



30.0 
= 14.8 
= 0.93 
= 138 



bot 



0.1 




50 100 150 50 100 150 

r (mm) r (mm) 

(a) Out-of-plane displacement {w£). (b) In-plane displacement (u t ). 



Figure 7 - Comparison of displacements from one-dimensional and two- 
dimensional finite element simulations for the model in Tablef2~\ 
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6 S UMMARY AND CONCL USIONS 



Example #2: Core Shear Stress 
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Example #2: Transverse Stress - Top of Core 
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(b) Transverse stress (a zz )- core top. 



Figure 8 - Comparison of stresses from one-dimensional and two-dimensional fi- 
nite element simulations for the model in Table^ 
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Example #2: Transverse Stress - Bottom of Core 
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Figure 9 - Comparison of stresses and displacements at the bottom of the core 
from one-dimensional and two-dimensional finite element simulations 
for the model in Table^ 
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6 Summary and Conclusions 

A detailed on-dimensional theory for sandwich panels with inserts has been derived. The approach 
follows that used by Thomsen (H. The models has been discretized using a finite element approach. 
The one-dimensional model produces results that are close to those of a two-dimensional 
axisysmmetric finite element model. Both models assume that the core is homogeneous, indicating 
that the one-dimensional model might be well suited for small deformations of sandwich specimens 
with foam cores. Further work is need to find nonlinear one-dimensional models of sandwich panels 
with inserts. 
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